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ABSTRACT 

Recent spectroscopic observations of transiting hot Jupiters have permitted the derivation of 
the thermal structure and molecular abundances of H2O, CO2, CO, and CH4 in these ex- 
treme atmospheres. Here, for the first time, we apply the technique of optimal estimation 
to determine the thermal structure and composition of an exoplanet by solving the inverse 
problem. The development of a suite of radiative transfer and retrieval tools for exoplanet 
atmospheres is described, building upon a retrieval algorithm which is extensively used in 
the study of our own solar system. First, we discuss the plausibility of detection of different 
molecules in the dayside atmosphere of HD 189733b and the best-fitting spectrum retrieved 
from all publicly available sets of secondary eclipse observations between 1 .45 and 24 /im. 
Additionally, we use contribution functions to assess the vertical sensitivity of the emission 
spectrum to temperatures and molecular composition. Over the altitudes probed by the con- 
tribution functions, the retrieved thermal structure shows an isothermal upper atmosphere 
overlying a deeper adiabatic layer (temperature decreasing with altitude), which is consis- 
tent with previously-reported dynamical and observational results. The formal uncertainties 
on retrieved parameters are estimated conservatively using an analysis of the cross-correlation 
functions and the degeneracy between different atmospheric properties. The formal solution of 
the inverse problem suggests that the uncertainties on retrieved parameters are larger than sug- 
gested in previous studies, and that the presence of CO and CH4 is only marginally supported 
by the available data. Nevertheless, by including as broad a wavelength range as possible in 
the retrieval, we demonstrate that available spectra of HD 189733b can constrain a family of 
potential solutions for the atmospheric structure. 

Key words: physical data and processes: radiative transfer - methods: analytical - planetary 
systems: individual: HD 189733b - planets and satellites: atmospheres 



1 INTRODUCTION 

The vertical structure of an exoplanetary atmosphere can be de- 
rived from transmission and emission spectroscopy when a planet 
transits its host star. The atmospheric properties of HD 189733b 
have been extensively investigated using transit spectroscopy ob- 
servations from both space and ground-based telescopes. Trans- 
mission spectra acquired during primary transits, where light from 
the host star is filtered through the upper layers of an exoplane- 
tary atmosphere, have been measured in the near-, mid- and far- 
infrared (NIR, MIR and FI R, respectively) by both the Hubble 
Space Telescope (H ST) JSwain et alj [2008: ISing et al. ' " 



and Spitzer teles cope dTinetti et alj 2007 J Beaulieu et al.1 



2009) 



2008 



IDesert et alj|2009h . Spitzer mid- and far-infrared (MIR and FIR) 
broadband photometry was used to deduce a high abundance of 
H2O in the terminator region of this planet jTinetti et al" 2007 al ; 
ISwain et al ] |2008| ; iBeaulieu et ail 120081 : IDesert et alj l2009h . In 
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addition, a strong absorption due to CH4 was reported in the 
near-infrared (NIR) using the NICMOS (Near Infrared Camera 
and Multi-O bject Spectrometer) spectrophotometry on the HST 
dSwain et al . 2008). However, the accuracy of transit spectroscopy 
and photometry has be en called into questio n (see the rev i ew by 
ISeager & Demindl201Ch . and most recently iGibson et alj J201lh 
claimed that instrumental systematics limit our ability to constrain 
atmospheric CH 4 from these HST/NICMOS data as well as C0 2 
and CO. A spectral feature at 4.5 /im can be explained by either 
of the carbon bearing molecules CO2 or CO, althou gh it is still 
unclear which due to the degeneracy of the solution ( IDesert et al.1 
2009; iFortnev et al.ll2010D . At visible wavelengths, strong lines of 
the alkali metal sodium were reported using observations with the 
high r esolution ground spe ctrograph of the Hobby-Eberly Tele- 
scope (Redfield et al. 2008), while featureless sp ectra detected by 
the HST/ACS jPont et alj2008l) and HST/STIS JSing et alj20"Tl}) 
are thought to be caused by thick atmospheric hazes. 

Observations of secondary transit spectra of HD 189733b by 
Spitzer broadband photometry, when emission spectra from the 
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planet's dayside is detected before and after the planet is eclipsed, 
have suggest ed that H2O is also abundantly present in the dayside 
atmosphere dCharbonneau et al.ll2008l : iDeming et ai1l200d) . form- 
ing the main features of the emission spectrum. The IrbO-rich at- 
mosphere was confir med via Spitzer IRS ( Infrared Spectrograph) 
dayside spectroscopy dGrillmair et ai1l2008» . #ST/NICMOS day- 
side observations in dicated a high abundance of CO2, CO, and CH4 
dSwain et alj|2009l. hereafter S09), supporting the conclusions of 
their primary transit spectrophotometry. With all measurements, 
each of these studies suggested that an atmospheric thermal pro- 
file lacking a thermal inversion (i.e. without a warm stratosphere) 
could best explain their observations. Finally, a decreasing temper- 
ature with altitude between 0.01-1.0 bar suggeste d the presence of a 
troposphere (S 09; IMadhusudhan & Seagerj|20 09. hereafter MS09). 
This work revisits the secondary transit spectroscopy of this hot 
Jupiter to assess the family of atmospheric structures consistent 
with these measurements. 

Various techniques have been used to investigate the at- 
mospheric constituents and structure from transit spectroscopy. 
Infrared dayside spectra during secondary eclipses measure the 
thermal emission of the dayside of the exoplanet, which de- 
pends mostly on the vertical temperature structure, while primary 
eclipse observations are more sensitive to the absolute te mpera- 
ture at the terminator due to its scale height dependence (Brown 
1200 lL iTinetti et al.l 12007b ;). The retrieval of atmospheric proper- 
ties from remot e sensing measuremen ts has been developed over 
many decades dGoodv & Yunel fl989) and has become a com- 
mon tool for the study of planetary atmospheres. Measured ra- 
diances (either photometry or spectroscopy) can be compared to 
synthetic spectra from radiative transfer modelling so as to de- 
termine the atmospheric structure and composition. Based on the 
Bayesian approach, an optimal estimation retrieval model solves 
an ill-constrained inversion problem using an iterative scheme to 
maximise the probability of solut ions to the av ailable data (the 
maximum a posteriori solution) (Rodgers 2000). Previous stud- 
ies of the t e rminat or and daysid e spec trum, such as those by 
ITinetti et al.l d2007ah . ISwain et alj d2008l) and S09, have utilised 
forward modelling to determine a best-fitted spectrum to the 
data, and have constrained the range of molecular abundances 
of transiting exoplanets by synthesising spectra based on theo- 
retical pressure-temperature (P-T) profiles. An alternative for- 
ward modelling approach used by a number o f studies ( MS09 
Sing et al.l 120081: IMadhusudhan & Seageri |2010|; IStevenson et al 
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used freely roaming P-T profiles and molecular abundances in 
a parameterised space, where calculated spectra were again com- 
pared with observations in terms of the goodness-of-fit, and their 
numerous runs were able to constrain P-T profiles and atmospheric 
compositions. Although these forward modelling techniques pro- 
vide valuable insights, their methods are based on line-by-line ra- 
diative transfer models, which are slow when large numbers of syn- 
thetic spectra are to be calculated over a wide wavelength range. 
Moreover, the huge degeneracies between the different model pa- 
rameters were not explored in detail. 

The NEME SIS optimal estimation retrie val algorithm 
I Irwin et all 1 2008) uses the correlated-fc technique dLacis & Oinaj 
199ll) in~its radiative transfer model, which rapidly integrates 
synthetic spectra using fc-distribution tables pre-calculated from 
line databases dGoodv & Yun3 ll989). The combination of the 
correlated-fc forward model, which is orders of magnitude faster 
than a standard line-by-line m odel, with an optimal estimation re- 
trieval scheme dRodgersl2 000) has been used to successfully inves- 



tigate planetary atmospheres in our own solar system. In this study 
we apply this rapid retrieval architecture to exoplanets for the first 
time, which allows us to formally address the uncertainties and de- 
generacies inherent in previous studies. Moreover, we use the co- 
variance matrices calculated via the inversion method to quantify 
the correlations between the derived P-T profile and atmospheric 
composition. Consequently, this study will show how the charac- 
teristics of these extreme exoplanet atmospheres can be deduced 
from spectroscopic measurements and highlight the limitations of 
the datasets available today. 

In this report, the modelling method used to retrieve the at- 
mospheric properties from the dayside spectroscopy and photom- 
etry of HD 189733b is described in Section 2. All available mea- 
surements used to constrain atmospheric properties are described 
in Section 3 and the best-fitting solution to these measurements is 
discussed in Section 4. Section 5 demonstrates the validity of the 
retrieval scheme, and quantifies the degeneracy of the solution. In 
Section 6 we discuss the implications of our results and arrive at 
our conclusions. 



2 MODELLING 

The atmospheric thermal structure and composition are de- 
rived using an op timal estimation retrieval algorithm, NEMESIS 
dlrwin et al]|2008h . Instead of simply minimising the deviation be- 
tween data and model by comparing the observations to millions 
of individually calculated synthetic spectra (the forward model), 
this algorithm instead solves the inverse problem using an iter- 
ative scheme. The retrieval problem is ill-posed (i.e. the number 
of unknowns is greater than the number of measurements) and ill- 
conditioned (i.e. small changes in the measurements can potentially 
have a disproportionately large effect on the solution) so there is no 
single unique solution. Instead, physically-realistic solutions and a 
range of u ncertainty must be defined for the family of possibilities. 
Following Rodgers (2000), we derive the optimal atmospheric state 
(x) via the inclusion of additional a priori constraints, in the form 
of an a priori state vector, x a , with error covariance matrix, S a . 
T hus successive iterations are used to minimise the cost function, 
4> dlrwin et al.l2008l ; lRodgersll2000l) : 

4> = (y - F(x)) T sr 1 (y - F(x)) + (x- x a ) T S- 1 (x - x a ) 

where F(x) is the forward model (i.e. the spectrum calculated for 
the trial atmosphere); x is the state vector (which may be a continu- 
ous atmospheric profile of temperature or composition, or a scaling 
factor of a known profile) and y is the measurement vector with 
measurement error covariance matrix, S e . The first term requires 
that the synthetic spectra must give a close fit to the measured spec- 
trum, th e second term r equires that the solution remains physically 
realistic. Rodgers (2000) shows that the optimal estimation solution 
for the state vector is given by: 

x = x a + (k t S7 1 K + S" 1 ) _1 K^S^y - Kx a ) 

where K is the matrix of functional derivatives (i.e. the rate 
of change of radiance with each model parameter). The op- 
timal state in the algorithm is accomplished after a moderate 
number of iterations (approximately 10-20), using a scheme 
based on the M arquardt-Levenberg principle (Levenberg 1 19441 ; 
lMarquardt|[l963l) . The algorithm calculates synthetic spectra using 
the correlated-fc approximation method, which has been shown to 
be fast, reliable, and sufficiently accurate compared to line-by-line 
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calculations. We compared spectra computed by line-by-line and 
correlated-fc models, and found that the mean difference between 
the two techniques is less than 5 per cent, which is considerably 
smaller than the range of uncertainties on available exoplanet mea- 
surements. 

In the fc-distribution method, the absorption spectrum over an 
interval is sorted in order of increasing absorption and the fraction 
of the interval with absorption less than a certain value k(g) is rep- 
resented in terms of the fraction of the interval g. Since k(g) is a 
smoothly varying function of g, it may be integrated with relatively 
few quadrature points N to determine the mean transmission over 
the interval as: 

JV 

T(m)^^e- fe ^A 5l , 

i=l 

where m is the amount of a molecule; and hi and Agi are the 
fc-coefficients and quadrature weights at each quadrature point 
dlrwin et ail 2008). The fc-coefficients are calculated from line data 
for a range of temperatures and pressures expected in the atmo- 
sphere in advance of the retrieval to enable rapid calculation of the 
transmission during each iteration. The mole cular line lists used in 
this study were taken from HITEMP2010 dRofhman et"ai]|2010l) 
for H2 O, the Carbon Dio xide Spectroscopic Databank (CDSD- 
1000) dTashkun et alj2003l) (also use d for HITEMP2010) for C0 2 , 
HITEMP1995 dRothman et aljl995l) for CO, and th e Spherical Top 
Data System (STDS) dWenger & Championlll998l) for CH 4 . The 
absorption coefficients for the extreme atmospheric temperatures 
found within exoplanetary atmospheres are continuously being im- 
proved, so the conclusions of this study are considerably limited by 
the quality of the available spectroscopic data. We have made use 
of the best available spectroscopic parameters at the time of writ- 
ing, and our line database will be updated as new sources of data 
become available. 

The combination of the optimal estimation retrieval scheme 
and the correlated-fc method enhances the efficiency of the retrieval 
process in terms of time and computational resources. Furthermore, 
NEMESIS calculates the matrix of the partial derivatives of radi- 
ances at each wavelength with respect to each retrieved variable, 
which are called the functional derivatives (or Jacobians), in 
order that the contribution of different atmospheric parameters at 
each wavelength can be easily interpreted by comparing the ele- 
ments of this matrix. 

Our a priori dayside atmosphere of HD 189733b extends 
from 10~ 9 to 10 bar to capture the full atmospheric range of poten- 
tial spectral contributions. For initial modelling we assumed that all 
species were well-mixed throughout the atmosphere and the molec- 
ular abundance was defined in terms of a single scaling parameter. 
This is because the retrieval of a continuous profile of composition 
would be under-constrained considering the small number of data 
points and the low spectral resolutions available, leading to non- 
physical oscillations in the retrieved profile. The a priori estimate 
for the abundance scaling parameter is assumed to have a large un- 
certainty so that retrieved values are not weighted by initial guesses 
since a simple scaling parameter already includes vertical smooth- 
ing. The a priori for temperature, however, is assumed to have a 
continuous profile and the assumed error on the a priori profile 
was adjusted to achieve the optimal balance between the quality of 
the fit to the measured data and the vertical smoothing. This will be 
further discussed in Section 5.3.1. 

As an important source of absorption in exoplanetary at- 
mospheres, collisional-induced absorption (CIA) between princi- 



ple gases are included in the model atmosphere. We consider 
the interactions between H2-H2 and H 2 -He and the coefficients 
are taken from iBorvsow et al.l (l989l). iBorvsow & Frommholdl 
dl989h.lBorvsow & Frornn^olddl99d).|Zheng & Borvsowldl995E 
IBorvsow et al. ( 19971) , and lBorvsowl d2002l) . The mole fractions of 
H2 and He are assumed to be related to the fractions of atomic H 
and H e, which are close to the typical solar value of 0.91 and 0.0887 
each dBurrows & Shard 19991) . The implications of this assumption 
will be tested in Section 5.5. 

Scattering by clouds and hazes has been reported from HST 
observations at visible wavelengths, yielding a featur e less trans- 
missi on spectru m of HD 189733 b dPont et alj [20081 ; ISing et al.l 
l201lh . Recently. Illeng et al.l J201 lh showed that the effect of scat- 
tering by clouds and hazes modifies the inferred temperature pro- 
file, and that an isothermal temperature structure above an adiabatic 
troposphere could also be caused by a cloud-top or haze layer. We 
consider, however, that the inclusion of scattering clouds and hazes 
in our retrieval model is beyond the scope of the present study for 
a number of reasons: (i) the scattering properties of possible clouds 
and hazes is insufficiently known; (ii) adding such particles would 
greatly enlarge an already large and poorly constrained parameter 
space and; (iii) modelling scattering processes will significantly in- 
crease the computation time. In summary, the addition of scattering 
would increase the complexity of the model and does not appear to 
be warranted by the dayside emission spectra studied here. We in- 
tend to assess the likely effects of scattering on our modelling of 
hotjupiters in a follow-up study. 



3 DATA 

To retrieve the P-T profile and compositional abundances, we 
used three measurement sets of the secondary eclipse for HD 
189733b, representing all of the available infrared dayside eclipse 
measurements at the present time, ranging over a wide wave- 
length range from 1.45 am to 24 um: (i) eighteen i/ST/NICMOS 
(Swai n et al.ll2009h chann els covering the rang e 1.45-2.5 pm; (ii) 
forty seven Spitzer IRS dGrillmair et ai]|2008b channels coverin 
the range 5-14.5 urn and one IRS photometry dDeming et 



covering 
alj2006t) 



channel at 16 um; and (iii) five Spitzer IRAC (Infrared Array 
Camera) and MIPS (Mult iband Imaging Photometer for Spitzer) 
dCharbonneau et al.l [2008) photometry channels at 3.6, 4.5, 5.8, 
8.0, and 24 um. The measurement errors on the three observa- 
tional datasets are directly obtained from the studies referenced 
above. To integrate the binned flux in each channel, we use (i) fil- 
ter widths taken from the literature for Spitzer broad band chan- 
nels dFazio et alj[2004l ; iRieke et alj|2004l ; IPeming et alj|2006l) : (ii) 
7J5T/NICMOS widths of -10 nm at 2 um; and (iii) Spitzer /IRS 
widths of — 100 nm at 8 (jm. The reference stellar spectrum for HD 
189733 is taken from the Kurucz grid modefl 

Despite the substantial efforts towards finding and reducing 
the errors from the data, uncertainties on the given datasets still re- 
main and are widely distributed over the wavelengths due to various 
error sources. First of all, techniques to decorrelate a transit light 
curve from the combi ned light of a plan etar y system are not cq n- 
sistent each other re.g. lSwain et al.ld2008l) vs. lGibson et al.ld201lh l. 
We will discuss the effects of underestimated errors on the dayside 
spectra in Section 4.1. Also, systematic errors inherent in the ob- 
servations themselves could be a strong error source. These two 



http://kurucz.harvard.edu/stars/HD 1 89733/ 



4 J. -M. Lee, L. N. Fletcher and P. G. J. Irwin 
Dayside emission spectrum of HD 189733b 



o Spitzer IRS (Grillmair et al. 2008) 
□ Spitzer IRAC/MIPS (Charbonneau et al. 2008) 
■ Spitzer IRS broadband (Deming et al. 2006) 
o HST NICMOS (Swain et al. 2009) 
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Figure 1. The best-fitting spectrum and atmospheric properties retrieved by 
the NEMESIS algorithm. In the top panel, the best-fitting dayside emission 
spectrum of HD 189733b is displayed as a black line. The green, purple 
with red, and green symbols ar e the measu red planet-star flux ratio from 
the Spitzer IRS spectroscopy jGrillmair e t al. 2008), the Spitz er broad- 
band photometry ( Charbonneau et al. 2008 c iDeming et al]|2O06l) . and the 
HST/NICMOS spectrophotometry JSwain et all2009h 7c7 Section 4). The 
retrieved best-fitting P-T profile in the bottom-left panel is shown in a black 
solid line with grey-coloured uncertainties due to the molecular degeneracy 
(cf. Section 5.3). In the bottom-right panel, retrieved best-fitting molecular 
mixing ratios for H2O, CO2, CO, and CH4 are demonstrated with differ- 
ent colours. Each abundance is assumed to be fixed with height in order to 
minimise the number of retrieval variables (cf. Section 5.4). 



4 BEST-FITTING DAYSIDE SPECTRUM OF HD 189733b 
4.1 Significance test (F-test) 

The molecular composition of the dayside atmosphere of HD 
189733b is generally accepted to be a mixture of H2O, CO2, CO, 
and CH4. These molecules have been discovered based on var- 
ious spectroscopic measurements, listed in Section 3. However, 
given the large uncertainties on the available data and the de- 
generacy between retrieved parameters, the existence of some of 
these species is called into question. We therefore set out to sta- 
tistically test whether the presence of each molecule is truly war- 
ranted by the data. Using a statistica l significance test known as 
an F-test dSnedecor & Cochran| [l989). we evaluated the change in 
confidence level between a range of simple and complex models. 
We successively added more complexity to the model (i.e. added 
molecules to the H2/He atmosphere to increase the number of de- 
grees of freedom) and evaluated the improvement in the goodness- 
of-fit parameter (x 2 ) for the best-fitting spectrum. The _F-test was 
then used to assess whether the inclusion of a particular molecule 
was required to fit the data. 

The F-test indicates that an atmospheric composition for HD 
189733b containing H2O and CO2 are highly plausible at the 
>99.98% confidence level. The presence of these two molecules 
is required to obtain a reasonable fit to the available data. In addi- 
tion, we find insignificant confidence levels (i.e. <95%) when we 
increase the complexity of the model by adding CO, CH4 or a com- 
bination of the two to the simple model containing H2, He, H2O 
and CO2 . In other words, CO and CH4 provide negligible enhance- 
ments to the fitting quality and are not required to fit the dayside 
emission measurements to within the stated error bars (although 

upper limits on these species can certainly be derived). 

Additionally, we consider the implications of iGibson et al.l 

(2011), who claimed larger uncertainties on t he HST/NICMO S 
data by re-processing its transmission spectrum ( Swain et al . 2008). 
If the same conclusions are applicable to the secondary eclipse 
emission spectrum, then we must similarly increase the measure- 
ment error on the NICMOS data by a factor of five. This leads to 
even smaller confidence levels on the more complex models, and 
the solutions are found to be even more degenerate. In this case, 
evaluating the F-test significance using a variety of models still 
suggests that both H2O and CO2 are required to reproduce the day- 
side emission spectra (>99.29% confidence), but makes the pres- 
ence of CO and CH4 even more uncertain. 

We therefore conclude that, irrespective of the uncertainties 
on the HST/NICMOS data, current datasets are unable to provide 
detections of CO and CH 4 on the dayside of HD 189733b with 
any sort of statistical certainty. Nevertheless, upper limits on the 
abundances of these molecules can be obtained (see Section 5), and 
all four molecules will be included in our subsequent study. 



facts may cause the planet-star flux ratio to be dependent on the 
data reduction process. Non concurrent observations for different 
wavelengths could cause large variations in the measured flux be- 
tween observations if there is any significant temporal variability 
in the atmosphere of HD 189733b. Therefore, all these sources of 
potential error could lead to substantial inconsistencies between the 
datasets taken from different studies, instruments, and observation 
times. Although we do not intend to revisit the decorrelation pro- 
cedures with this analysis, we aim to show the limitations of the 
spectra available at the present epoch under the assumption that the 
measurements are an accurate representation of the mean planetary 
flux. 



4.2 Best-fitting solution 

Using the NEMESIS algorithm, we retrieve the best-fitting day- 
side spectrum of HD 189733b, incorporating both the Spitzer and 
HST observations as stated in Section 3. Fig.[T]demonstrates the 
best-fitting spectrum to these measurements, and the retrieved at- 
mospheric P-T profile and molecular abundances for H2O, CO2, 
CO, and CH4. Each panel will be described below. 

Fig. [2] shows the contributions from the four main gases in- 
cluded in our model. This figure also shows the wavelength ranges 
where the molecular contributions are distributed by co-plotting the 
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Figure 2. Fitted dayside emission spectra of HD 1 89733b. All available 
a black line in all figures. We also show calculated spectra with various 
best-fitting spectrum. For H2O, molecular abundances are varied by 0.1 
0.01, 0.05, 20.0 and 100.0 times for the other molecules. 



measurements are shown in grey symbols. The best-fitting spectrum is displayed as 
molecular abundances to understand the contributions of different molecules to the 
, 0.2, 5.0 and 10.0 times the abundance leading to the best- fitting spectrum, and by 
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computed synthetic spectra with high and low abundances for each 
molecule. Spectral features of H2O and CO2 affect the spectrum 
at all wavelengths. CO has absorption features at 1.6 —1.7, 2.3-2.5, 
and 5.0-5.5 /im, while CH4 features can be seen at 1.7, 2.2-2.5, 
3-4, and 7-9 pm. A striking feature of the best-fitting spectrum 
is a deep IR absorption by CO2 at 9-24 ^m, which is rather dif- 
feren t from the fitting by other studies (MS09;|Fortnev & Marled 
120071 : iGrillmair et alTl2008l ; ICharbonneau et alj|2008h . who con- 
cluded that the features at the longer wavelengths of IRS are caused 
mainly by H2O, of which absorption is also dominant at the shorter 
wavelengths of 5-9 pm. Although the low CO2 hypothesis seemed 
to give a good reproduction of the long-wave Spitzer data, it fails 
to reproduce the shorter wavelengths covered by HST. As seen 
here, one of the benefits of our approach is that the use of a broad 
range of IR wavelengths allows us to break some of the degen- 
eracies inherent in modelling a small number of data points (cf. 
Section 5.4). 

In this study, the high CO2 abundance causes a sharp drop at 
9 /im and then a flat and featureless spectrum between 9-24 pm. 
H2O still contributes to the spectrum at 5-9 and >20 pm as 
well as a small feature from a vibrational transitio n at 6-6.5 /-tin, 
which is one of the H2O features highlighted by IGrillmair et al.l 
(2008). In particular, the spectrum fits the planet-star flux ratio for 
the Spitzer IRAC 3.6 pm channel, a feature which has been previ- 
ously explained by the large 3.25 pm emission of CH4 in non-local 
thermodynamic eq uilibrium (NLTE) conditions dSwain et alj20ld : 
iThatte et alj|2010l) . However, any additional NLTE calculation is 
not required in our model to produce the data at 3.6 ^m adequately. 
For the HST/NICMOS channels between 1.45-2.5 ^m, we fit to 
the measurements using a high amount of CO2 as suggested in 
MS09 (~7x 10~ 4 ) rather than in S09 whose model suggested low 
CO2 mixing rati o (10~ 7 -10~ 6 ). This discrepancy is not yet fully 
resolved, but, as Shabram e t alj d201 lh stated, this may come from 
the difference in the forward models used for each study. In our 
spectrum, the channels at 1.584, 1.869, 2.159-2.216 ^m, however, 
are underestimated by the high CO2 abundance, leading to a de- 
crease in the quality of the fit to the data. The sharp edge present in 
our synthetic spectrum at 9 pm, w here a sudden CO2 line weaken- 
ing is shown dTashkun et afl |2003). is caused by a strong CO2 ab- 
sorption. Its absence from the measured exoplanet spectrum is yet 
to be understood but may be related to the limited number of gases 
(H2O, CO2, CO, and CH4) and the absence of absorbing aerosols 
considered in our modelling. 



5 RETRIEVAL OF ATMOSPHERIC PROPERTIES 

On the basis of the best-fitting spectrum presented in the previous 
section, we here validate the retrieval method by analysing the tem- 
perature contribution functions and the functional derivatives for 
the molecular abundances. These sensitivity analyses then allow 
plausible interpretations for temperature and composition structure 
of the dayside of HD 189733b. 

5.1 Contribution functions 

The radiance measured in each channel originates from different 
pressure levels within the atmosphere. This is generally described 
using a contribution function, which is the product of the Planck 
function at the local temperature in the atmosphere and the trans- 
mission weighting function, which describes the rate of change of 
atmospheric transmission with height. The contribution functions 



are dependent on the best-fitting profiles for the atmospheric struc- 
ture and abundances. Thus the calculated contribution functions in- 
dicate the pressure levels at which thermal emission from the at- 
mosphere contributes most to the radiance observed in each chan- 
nel. Fig. [3] shows the contribution functions for all channels of the 
Spitzer IRS, IRAC, and MIPS, and HST/NICMOS. 

The six Spitzer broadband photometry channels in the range 
3.6-24 pm have broadly-distributed contribution functions whose 
peak pressures range from 2 to 300 mbar. The contribution func- 
tions for the IRAC channels (3.6, 4.5, 5.8, and 8 pm) and MIPS (24 
pm) are located in the deeper atmosphere and provide strong con- 
straints for the temperature between 30-300 mbar. The contribution 
function of the IRAC 4.5 pm channel has a second peak at high al- 
titude (3-5 mbar), being close to the peak of the IRS 16 pm channel 
at 2 mbar. Hence, the temperature at pressures as low as 2 mbar can 
be retrieved from the IRAC channels, but only if the temperature 
of the deep atmosphere is well constrained from the other mea- 
surements. The 47 Spitzer IRS spectroscopy channels between 5 
and 14.5 pm have closely-spaced overlapping contribution func- 
tions, with peaks ranging from 2 to 150 mbar. The radiance from 
the lower atmosphere (~100 mbar) contributes to the channels in 
the range 5-9 /im, whereas radiance from the upper atmosphere 
(>10 mbar) contributes longward of 9 pm. Unlike the MIR and 
FIR channels on Spitzer, the HST/NICMOS NIR 1.45-2.5 ^m 
channels can measure the emission from the deeper atmosphere at 
20-600 mbar and their contribution functions are partly distributed 
over the pressure levels that are not covered by the Spitzer mea- 
surements. 

In summary, the contribution functions show that the P-T pro- 
file over a broad range of pressures from 2 to 600 mbar can be con- 
strained by retrieving from the Spitzer and HST measurements 
together. On the other hand, studies that focus on a single narrow 
wavelength range are sensitive to only a narrow altitude range, with 
potential degeneracies between temperatures at high and low alti- 
tudes. Thus, by considering all available datasets simultaneously, 
we can provide stronger constraints on the retrieved thermal prop- 
erties than in previous studies. 

5.2 Functional derivatives 

The functional derivatives are defined to be the partial derivatives of 
the radiance (or any spectral output in the forward model) with re- 
spect to any given atmospheric parameters. By calculating the func- 
tional derivatives for molecular abundances, we can understand 
which measurements are sensitive to the molecular abundance and 
at which pressure levels. Fig.|4]shows the functional derivatives for 
the 4 molecules considered in this study. For ease of comparison, 
these are normalised to the peak of the functional derivatives for 
each measurement set. Using Fig. [4] the sensitivity can be inter- 
preted in two directions: Each row indicates which pressure levels 
show high sensitivity to a given molecule in each channel, and each 
column indicates how molecular sensitivity is distributed through 
pressure levels in a given set of measurements. In all cases, the 
derivatives are negative because an increasing abundance increases 
the atmospheric absorption (in the absence of a thermal inversion 
on this planet) and decreases the ratio of the disk-averaged plane- 
tary flux with respect to the host star. 

High sensitivity to the abundance of H2O is seen in all 
channels, with the exception of the 1.93-2.16 pm channels of 
HST/NICMOS, 4.5 channel of Spitzer IRAC, 9-14.5 pm 
and 16 channels of the Spitzer IRS, in which regions the mod- 
elled radiance is dominated by CO2 absorption, as explained below. 
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Figure 3. Contribution functions for the Spitzer broadband photometry (left), IRS spectroscopy (middle), and the H ST/NICMOS spectrophotometry (right) 
channels. For the Spitzer photometry channels, each line pattern means MIPS 24 fim (solid), IRS 16 fim (solid-black), IRAC 8.0 /im (dotted), 5.8 fim 
(dashed), 4.5 pm (dot-dashed), and 3.6 pm (triple dot-dashed). For the Spitzer IRS and the H 5T/NICMOS channels, the brighter colours denote the 
channels at the shorter wavelengths. For all cases, emission from the lower atmosphere tends to dominate the shorter wavelength channels. 




Figure 4. Normalised functional derivatives of molecules in Spitzer and HST channels. The applied colours are the same as those described in Fig. [3] Each 
row shows the vertical sensitivity of radiance with respect to the abundances H2O, CO2, CO, and CH4 (from top to bottom). This figure shows the pressure 
levels at which molecular abundance can be retrieved (row direction) from which channel (column direction). 
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Figure 5. (a) Retrieved P-T profiles of the day side of HD 1 89733b from a 
range of diverse a priori profiles. The different lines show: a priori P-T 
profiles (dashed), retrieved P-T profiles (solid) and their errors (dot-dashed 
and dotted). This plot shows the pressure range over which the P-T profile 
is retrievable from these measurements, (b) Retrieved P-T profiles from 
the same a priori profile shape, but offset from each other with baseline 
temperatures at 1600 K, 2000 K, and 2400 K at 10 bar. 

The H2O functional derivatives all peak in the ~ 100-500 mbar re- 
gion, showing that the measurements can only constrain the H2O 
abundance in this altitude range. In contrast, the CO2 functional 
derivatives are divided into two separate pressure levels (peaks be- 
tween ~0. 1 bar and 0. 1-1 mbar) as are clearly shown in the second 
row of Fig.H The ffST/NICMOS channels between 1.584-1.698 
and 1.929-2.216 /im are only sensitive to the CO2 abundance at 
altitudes below 100 mbar whereas the Spitzer broadband photom- 
etry and IRS spectroscopy are sensitive to the CO2 in the 0.1-1 
mbar region as well. Despite this sensitivity to a range of altitudes 
for CO2, a combination of these channels, however, may not deter- 
mine CO2 abundance at both levels due to the small sensitivity in 
the upper atmosphere (cf. Section 5.4). CO can only be detected in 
three HST/NICMOS channels in the range 2.33-2.45 pm at 100 
mbar, and these channels have been extensively used to constrain 
the CO mixing ratio (S09 and MS09). The sensitivities for CH 4 are 
so marginal, as seen in the fourth row, that constraining its abun- 
dance may not be possible from the given measurements. 

In summary, the functional derivatives indicate the altitudes 
and channels showing sensitivity to each particular molecule. The 
sensitivities to the molecules are mostly clustered in the lower at- 
mosphere (~ 100 mbar), and, in particular, CO2 shows an additional 
peak at ~ 1 mbar. Therefore the molecular abundances can only be 
constrained in the deep atmosphere (troposphere), with less sensi- 
tivity to the upper atmosphere. Finally, these functional derivatives 
show that the spectrum is sensitive to the H2O abundance over a 
broad range in most of the measured channels, implying an inherent 
degeneracy between temperatures and H2O in these observations. 

5.3 Retrieval of P-T profile 

5.3.1 The best-fitting P-T profile 

Previously P-T profiles of HD 189733b have been estimated based 
on theoretical models or parametric retrievals to generate many 



thousands of model spectra to compare with observations. How- 
ever, this forward modelling approach does not explicitly solve the 
inverse problem during the constraining process and thus it is un- 
clear if the solutions are biased more towards theoretical expecta- 
tions than being driven by the measurements themselves. For this 
reason, we derive temperatures using several different a priori 
profiles to show that the retrieved temperature converges to a re- 
producible profile in the altitude range covered by the contribu- 
tion functions. Because a biased P-T profile would vary with the a 
priori assumptions, we look for an appropriate a priori error and 
its vertical shape in order that the same P-T profile can be retrieve 
irrespective of the shape of the a priori profile. Fig.[5ja) presents 
the retrieved P-T profile and its error for a range of selected tem- 
perature a priori. For all cases, the P-T profiles share a common 
shape between 0.1 mbar and 1 bar, demonstrating the validity of 
the temperature retrieval over this range. It is shown that even with 
the simplest possible assumption such as an isothermal temperature 
throughout the atmosphere [blue line in Fig. 0a)], the measure- 
ments still produce a similar thermal profile as the other retrievals. 
As a further test, we take an a priori structure from the retrieved 
profile in the previous step, offset it by a temperature 400 K, and 
repeat the temperature retrieval again. Fig.[5jb) shows that the P-T 
profile is still sufficiently constrained by the measurements, even if 
there are large shifts at levels not covered by the contribution func- 
tions, where the solutions relax back to their different a prioris. 
We conclude here that the available dayside emission spectra can 
constrain an unbiased atmospheric structure over the levels probed 
by the contribution functions. 

As a result, we find that the temperature decreases adia- 
batically from 1900 K at 600 mbar to ~1200 K at 100 mbar, 
then becomes isothermal up to the upper atmosphere (~1 mbar). 
These adiabatic and isothermal layers in the thermal structure 
are dominant features of heat transfer by convection and ra- 
diative cooling, respectively. In comparison S09 claimed a de- 
creasing temperature layer between 0.01 and 1 bar to model the 
H ST/N1CMOS measurements, making an adiabatic layer ~10 
times thicker than our estimation, and theoreti cal models also con- 
sidered adiabatic layers in the tropospher e dFortnev et al.l 1200^ : 
iBurrows et all 120081 : Ishowman et alj|2008t) . The isothermal tem- 
perature (~1100 K) of the dayside hemisphere at pressure levels 
above the troposphere (100 mbar) is possibly maintained by ef- 
ficient energy re-distribution throughout the whole planet system 
as suggested by the observations in the Spit zer IRAC channels 
dKnutsonet al. 2007; Charbonneau et al. 2008). Vari ous cir culation 
models dShowman et alj|2008l ; iDobbs-Dixon & Linll20 08) predict 
that the heat transports between the dayside and nightside lead to 
the isothermal structure at mid and low pressure (<100 mbar) in 
the dayside atmosphere. 

In summary, the broad wavelength range of available measure- 
ments has provided a strong constraint on the vertical P-T pro- 
file (an adiabatic troposphere and isothermal stratosphere) without 
excessive sensitivity to the a priori assumptions, and without re- 
liance on theoretical modelling or parameterised profiles. 

5.3.2 Temperature degeneracy 

In general, the best-fitting P-T profile from the retrieval is a non- 
unique solution. This means that there is potentially a large family 
of solutions for the P-T profile that could fit the measured spec- 
trum equally well within the same error range. In many retrievals, 
there may also exist cross-correlation with other fitted quantities, 
indicated by covariances that are larger than zero. The covariance 



HD 189733b atmospheric structure retrievals from transit spectroscopy 9 



(a) cross-correlation 



(b) T vs. all molecules 





J | 

i ; 

I ; 
1 : . 

^ : 

i ■■ 




l : 
f\\ 




V > \ 
: \\ 


- T(P)-H 2 


•\\ \ 


T(P)-C0 2 


'■■ 1 \ 


i.i \ 


T(P)-CO 


ii 

'<'''■ / 


T(P)-CH 4 


Hi / 




J . / 




S ; / 
'■'/ 




-0.5 



0.0 
c(i,j) 
(c) T vs. H 2 



0.5 



1000 2000 
Temperature (K) 

(d) T vs. C0 2 




1 000 2000 
Temperature (K) 

(e) T vs. CO 



1000 2000 
Temperature (K) 

(f) T vs. CH 4 




1 000 2000 
Temperature (K) 



1000 2000 
Temperature (K) 



Figure 6. (a) The vertical structure of the cross-correlation functions. The 
high correlation with H2O and CO2 dominates the degeneracy of temper- 
ature at most pressures, (b-f) The retrieved P-T profiles of the dayside 
HD 189733b with various mixing ratios for each molecule are presented 
with A\ 2 /N ranges. Each colour demonstrates A\ 2 /N <0.5 (red), < 1.0 
(green), and <2.0 (blue), respectively. As a result, the large uncertainties 
of P-T profile are broadly distributed between the deep and middle atmo- 
sphere. On the other hand, narrow uncertainty at 2 mbar implies an isother- 
mal structure in the upper atmosphere of the dayside of HD 189733b. 



matrix is advantageous because the matrix is computed as a part of 
the retrieval process and can be used in assessing the degeneracies 
between variables at the end of retrieval. 

For the dayside spectrum of HD 189733b, the degener- 
acy between the molecular abundances and the P-T profile is 
thought to be significant at some pressure levels and bandpasses 
(S09 and MS09), however, a detailed analysis of these correla- 
tions has not previously been presented. Here, we examined the 
cross-correlation functions c(i,j) (i.e. the off-diagonal elements 
of the correlation matrix) that determine the degree of degener- 
acy between elements i and j of the state vector. By definition, 
|c( i, j ) |=1.0 represents a perfect correlation between variable i and 
j, and we consider that | c|=0.5 is a limitation for quasi-independent 
retrievals of different variables. The profiles in Fig. |6ja) show the 
vertical structure of cross-correlation functions between the molec- 
ular abundances and the temperature at different pressure levels. 
It can be seen that the P-T profile is significantly correlated with 
molecular abundances at some levels, being most correlated with 
the H2O abundance between 200-400 mbar and with CO2 for alti- 
tudes above the 30 mbar pressure level. 

To determine the degeneracy of the temperature profile from 
Fig-UJa), we performed multiple retrievals where individual molec- 
ular abundances were fixed at a particular value and all other vari- 
ables retrieved. We then assessed the weighted mean squared er- 
ror, the x 2 over the number of measurements (N) of the solu- 
tion as a function of the fixed molecular abundance, compared 
with the minimum \ 2 /N (Xrnin/N) of the best-fitting retrieval, 
and determined the range of temperatures for which the change 
A. X 2 /N = x 2 /N ~ Xmin/N was less than 0.5, 1.0, and 2.0, re- 
spectively. 

Fig. |6lb) shows the ranges in the retrieved P-T profile for 
different values of A\ 2 /N due to the molecular degeneracy. At 
each pressure level, allowed temperature ranges are determined by 
varying the molecular abundance with which temperature is most 
strongly correlated. The temperature uncertainties for A\ 2 /N < 
1.0 calculated in this way are 540 K at 600 mbar, 420 K at 100 
mbar, and 100 K at 2 mbar. As described on the basis of cross- 
correlation, the temperature uncertainties at altitudes below 10 
mbar are caused mainly due to the degeneracy with the H2O abun- 
dance and those at altitudes above 10 mbar are degenerate with 
the CO2 abundance. On the other hand, the degeneracies between 
temperature and the abundances for CO and CH4 are small. There- 
fore it can be seen that although some uncertainty in the deeper 
P-T profile exists due to molecular degeneracy, the small degree 
of degeneracy at lower pressures supports our conclusion that HD 
189733b has an isothermal structure at the upper atmosphere. 

5.4 Retrieval of molecular abundances 

Despite the molecular degeneracies described above, we can use 
the dayside emission spectra to place crude estimates on the abun- 
dances in addition to the vertical temperature structure. Each gas 
is assumed to be well-mixed with a constant mole fraction with 
height. However, we know that the functional derivatives are only 
sensitive to a limited pressure range and thus our retrieved abun- 
dances represent mean concentrations only at the pressure levels 
which show the highest sensitivity (Fig. |4j. The best-fitting molec- 
ular abundances are determined from the retrieval accompanying 
the best-fitting P-T profile given in Section 5.3 and we find that 
those are 6.7 x 10~ 4 (H2O), 2.0x 10" 3 (CO2), 2.9 x 10~ 3 (CO), and 
1.9x 10 -7 (CH4), respectively. Due to its tiny contribution to the 
spectrum, the best estimation of CH4 should be considered as an 
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Table 1. Estimated and retrieved mixing ratios and C/O in the dayside of HD 189733b. 





H 2 (lO" 4 ) 


C0 2 (lO" 4 ) 


CO(10 -4 ) CH 4 (10- 4 ) 


C/O 


Data Source 


Line et al. 2010 


6-13 
~6.36 


0.0047-0.016 
0.004-~0.1 


2-9 0.0026-6.758 
-8.4 —0.4 




Thermochemistry 
Photochemistry 


Swain et al. 2009 


0.1-1 


0.001-0.01 


1-3 <0.001 


0.5-1 


.ffST/NICMOS 


Madhusudhan & Seager 2009 


0.01-1000 
0.1-10 
~1 


0.007-0.7 
~7 


<100 
<0.02 

2-200 <0.06 


0.007-1 
0.5-1 


Spitzer IRS 
Spitzer photometry 
HST/NICMOS 




6.7 


20 


Best-fitting value 
29 0.0019 


0.65 


All measurements 


This study 


0.9-50 
0.3-100 
0.03-400 


3-150 
1.5-300 
<2000 


Possible fit range 
<0.4 

<1 

<3 


0.45-1 
0.30-1 
0.15-1 


A X 2 /N < 0.5 
A X 2 /N < 1.0 
Ax 2 /N < 2.0 



upper limit only, and the CO abundance is poorly constrained by 
the available data (see below). Additionally, the C/O ratio of HD 
189733b (0.65) is found to be close to the solar value (—0.6). 

By looking at the covariance of our retrieved solution we find 
that the cross-correlation between the different molecular abun- 
dances is rather small (c < 0.1) for all the molecular combinations. 
However, there remains some degeneracy due to the main correla- 
tion with temperature. To quantify this we calculated A\ 2 /N for 
each gas by changing the abundance of a molecule and retrieving all 
other parameters including temperature. The variation of A\ 2 /N 
with abundance for all four gases is shown in Fig.[7] We find that the 
following abundance ranges have A x 2 /N < 0.5: (9-500)xl0~ 5 
for H 2 0, (3-150)x 10~ 4 for C0 2 , and <4x 10~ 5 for CH 4 , respec- 
tively. The error ranges for the abundance of CO were so broad that 
we cannot provide meaningful estimates from the dayside spectra 
alone, hence they are absent from Tab.[T] The retrieved abundances 
of the other species for A% 2 /N < 1.0 and 2.0 are also shown in 
Tab. Q] and are compare d with the exp ectations of a thermo- and 
photo-chemistry model Jline et aT]|2010l) . and with previous for- 
ward model studies (S09 and MS09). As we found using the F-test 
in Section 4.1, the data only constrain the upper bound of the mix- 
ing ratio for CH4 and do not offer useful constraints on the abun- 
dance of CO because the disk-averaged planet-star flux ratio shows 
only a small sensitivity to their abundances in Fig.|4]and their con- 
tributions to the emission spectrum appear in only few bandpasses. 
Consequently, variation of these two gases leads to a tiny change in 
A\ 2 /N where the abundances are lower than the best-fitting esti- 
mation. We conclude that it is not possible to constrain the abun- 
dances of CO and CH4 based on these datasets alone as previous 
studies have attempted. 

All channels of Spitzer and HST are capable of constrain- 
ing the H 2 mixing ratio as presented in the functional derivatives, 
and the constrained H 2 is consistent with other estimations. For 
CO, the strongest constraints are taken from the HSTTNICMOS 
data (MS09) and we find that the dayside emission spectra can 
offer no meaningful constraints on the abundance of this species. 
A relatively small abundance of CH4 has been estimated by pre- 
vious studies in Table 1, except for the large abundance (10 -2 ) 
derived from the Spitzer IRS spectroscopy at 7.6 fim (MS09), 
which is not consistent with this study. The derived abundance of 
C0 2 (3-150xl 0~ 4 ) is 2-3 orders of magnitude larger than pre- 
vious studies bv lLineetal.1 fcoidh (-10 -7 — 10 -5 ), S09 (10 -7 - 
10" 6 ), and MS09 (using Spitzer photometry, 7-700xl0 -7 ) for 




Volume Mixing Ratio 

Figure 7. The degeneracy ranges of the molecular mixing ratios for H 2 0, 
C0 2 , CO, and CH4. Each line shows resultant \ 2 /N with respect to given 
abundances. The A\ 2 /N constrains the uncertainties of the molecular 
abundances which are distributed around the best-fitting abundance of H 2 
(6.7xl0~ 4 ), C0 2 (2.0xK)- 3 ), CO (2.9xl0~ 3 ), and CH 4 (1.9x 10" 7 ), 
respectively. Lower bounds of CO and CH4 uncertainties are unconstrained 
because of their small contribution to the spectrum. 



all A X 2 /N ranges and thus our retrieved C0 2 abundance is only 
consistent with the ifST/NICMOS estimation by MS09. Indeed, 
it is the if ST/NICMOS channels that constrain our retrieved C0 2 
abundance to be larger than previous studies (— 10~ 3 ), whereas the 
Spitzer/IRAC channels at 4.5 and 16 /im could be fitted with ei- 
ther low or high abundances of this molecule. 

To test the sensitivity to C0 2 in the upper atmosphere, we 
compared the best-fitting spectra to test cases with a range of upper 
atmospheric abundances of this molecule (i.e. allowing the abun- 
dance to vary with altitude). The C0 2 mixing ratio was fixed at a 
certain value up to a pre-defined pressure level, and then declined 
linearly into the upper atmosphere with a slope determined by a 
fractional scale height (the ratio of the scale height of the gas to 
the scale height of the atmosphere). This allows us to demonstrate 
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Figure 8. Synthetic emission spectra of the dayside of HD 189733b with 
different CO2 abundances. The three lines show the spectrum with a well- 
mixed abundance of 2. Ox 10 -3 (black); well-mixed at 1 X 10 — 7 (cyan) and 
the case with CO2 varying with altitude from 2.0xl0~ 3 at 100 mbar to 
1.4x 10 -7 at 0.1 mbar (orange). The latter case is almost indistinguishable 
from the first scenario. 



the effect of different mixing ratios with height, but avoiding the 
need to introduce a complete vertical profile of CO2 (and hence a 
large number of additional parameters in our state vector). In Fig.[8] 
the retrieval with the new CO2 vertical profile that decreases from 
2.0x 10~ 3 at 100 mbar to 1.4 x 10~ 7 at 0.1 mbar, produces a very 
similar spectrum to the well-mixed high-C02 case, leading to an 
insignificant improvement in terms of the X*/N. This is because 
the CO2 sensitivity at the Spitzer bandpasses is too low to ad- 
equately constrain its abundance in the upper atmosphere (0.1-1 
mbar), whereas the strong constraints of the H STfNICMOS and 
Spitzer IRS channels between 9-16 /im allow us to constrain the 
abundance in the lower atmosphere (~100 mbar). In contrast, a 
small abundance of CO2 (10~ 7 ) with a well-mixed vertical profile 
is unable to reproduce the measured spectrum. Our results suggest 
that a large amount of CO2 (~10~ 3 ) exists at ~100 mbar but that 
the abundance at 0.1-1 mbar is unconstrained by current measure- 
ments. 



5.5 Additional Sources of Degeneracy 

In the sections so far, we have derived conclusions about the at- 
mospheric temperatures and composition based on a series of un- 
derlying assumptions, which have not been quantified by previous 



authors. Here we perform additional degeneracy tests to quantify 
the uncertainties on temperature and composition caused by the un- 
certainty of the amount of H2 and He and the radius ratio between 
planet and star (R p /R*). The abundances of H2 and He are gen- 
erally considered to be more or less the same as in the parent star. 
However, the helium abundance in the giant planets of our own so- 
lar system deviates considerably from solar, and this could provide 
an additional source of error in the retrieval of the vertical temper- 
ature structure. However, we find that the change in mole fraction 
of H2 (assumed to be l-q[He]) between 0.8-1.0 gives a negligible 
effect on the retrievals of temperature and composition from the 
dayside spectra at all pressure levels, suggesting that the emission 
spectrum of HD 189733b (which is dominated by the features of 
H2O) is not particularly sensitive to variations of H2 and He com- 
pared to the host star. 

On the other hand, variations of the planet-star radius ra- 
tio, which determines the bottom level of the model atmosphere, 
has a significant effect on the flux ratio. If we consider the effect 
of varying the ra dius ratio over the range 0.147<7? p /i?*<0.160 
dAgol et alj2oToh . then we find that the best-fitting molecular abun- 
dances vary over the ranges: H2O (6.5-6.7 x 10~ 4 ), CO2 (2.0- 
2.2x 10 -3 ), CO (2.5-3.3 x 10 -3 ), and CH 4 (l-10x 10~ 7 ), respec- 
tively (Fig.[9]l. The figure indicates that the temperature profile over 
the pressure range becomes ~100 K cooler when R p /R* is en- 
hanced by 0.013 (~10 per cent increase in the radius of the planet). 
As an increase in R p /R* enhances the flux ratio across all wave- 
lengths, the model tends to select a correspondingly low tempera- 
ture profile that decreases the flux ratio. 

Finally, as we described in Section 3, we have omitted the 
scattering and absorption processes due to the present of clouds 
and hazes in this hot Jupiter atmosphere, choosing to focus on a 
clear atmosphere for this initial exercise. Given the present qual- 
ity of the dayside emission spectra, the more complex model with 
scattering hazes is not needed to reproduce the data (although hazes 
are clearly needed to reproduce transmission spectra of the termi- 
nator region), and would lead to a significant rise in the number of 
model parameters. However, the omission of aerosols could lead 
to significant changes to our derived temperatures and molecular 
abundances, and our future work will aim to quantify these uncer- 
tainties to further explore the atmospheres of hot Jupiters. 



6 DISCUSSION AND CONCLUSIONS 

The retrieved thermal structure is a useful tool for further under- 
standing of dynamic processes in the atmosphere. As a descrip- 
tion of atmospheric stability, the lapse rate (r = dT/dz, the rate 
of change of temperature with respect to the height) is shown in 
Fig. [lO] An increase in the lapse rate indicates that there is an adi- 
abatic layer beneath the ~100 mbar level (~1400 km above from 
the level of 10 bar), and that the lapse rate tends to zero for pressure 
levels above ~30 mbar (~2100 km) where the temperature profile 
tends to an isotherm. Moreover, the measured lapse rate is similar 
to the dry adiabatic lapse rate (DALR) below the 300 mbar pres- 
sure level, and is overlain by a sub-adiabatic layer extending up to 
30 mbar, similar to the atmospheric thermal structures seen in the 
upper tropospheres of the giants planets of our solar system. 

The NEMESIS algorithm uses two novel tools for an efficient 
retrieval of the atmospheric properties: a non-linear optimal esti- 
mation scheme and a fast forward model taking advantage of the 
correlated-fc approximation and fc-distribution tables. This rapid re- 
trieval method allows us to infer general atmospheric properties 
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Figure 9. Uncertainty of composition (left) and temperature (right) being 
induced by the variation of planet-star radius ratio (R p /R t ). Due to ~10 
per cent increment in R v /R t , the temperature profile shifts towards lower 
~100 K whereas the fitting quality, A\ 2 /N increases by 0.5. 
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Figure 10. Disk-averaged dayside lapse rate of HD 1 89733b (right) from 
the derivative of the P-T profile (left). The given pressures correspond to 
the levels covered by the contribution functions. We also display the dry 
adiabatic lapse rate (DALR) in grey, assumed to be a constant over this 
pressure range. 



from small datasets. The optimal estimation retrievals permit an 
extension of previous results and the formal quantification of er- 
rors and uncertainties. The functional derivatives, which are ana- 
lytically computed, enable us to trace the vertical contribution of 
the different input constituents in each channel, permitting us to 
move beyond the single best-fitting approach and show the rela- 
tionship between the features of the data and physical values. We 
have shown that the current set of observations of the dayside spec- 
trum of HD 189733b are enough to constrain the thermal structure 
at some pressure levels and the mixing ratios for H2O and CO2, al- 
beit with large uncertainties due to molecular degeneracies and the 
limited number of spectral measurements. On the other hand, we 
do not find statistically-significant evidence for CO and CH4 in the 
dayside emission spectra alone, and the present observational set 
of HD 189733b only allows us to place upper limits on the abun- 



dances of these species. Furthermore, the ability to calculate the 
cross-correlation function allows us to assess the degeneracies be- 
tween the various modelled parameters in our state vector. 

With the retrieval method presented here, we derive three ma- 
jor findings from the measurements of the dayside of HD 189733b. 
Firstly, we have produced the first retrieval (i.e. solution to the in- 
verse problem) of the vertical temperature structure of the dayside 
atmosphere of HD 189733b. The retrieved profile shows that the 
constraints in the NIR and MIR lead to an adiabatically decreasing 
temperature between 0.1-1 bar. In addition, strong constraints from 
the Spitzer IRS measurements in the 9-16 /im range suggest that 
the temperature structure in the mid and upper atmosphere (1-100 
mbar) seems to be isothermal with a uniform temperature (-~1100 
K). One explanation for the isothermal layer is that super-rotating 
jets blowing from the night- to dayside play an critical role of an 
efficient energy re-distribution over the whole planet and may be 
responsibl e for maintaining the high temperature in the upper at- 
mosphere l lKnutsonet al. 2007; Showman et alj2008h . 

Secondly, the functional derivatives for the molecules show 
that the measurements at NIR (1.6-1.7 and 1.9-2.2 /im), MIR (4.5 
/im), and FIR (9-16 fim) spectral ranges constrain the CO2 abun- 
dance at ~100 mbar but not in the upper atmosphere. Its abun- 
dance in the lower atmosphere is responsible for the spectral flat- 
ness at 9-16 /im, where previous studies assumed H2O to be the 
dominant constituent. We determine the vertical sensitivity of CO2 
abundance by comparing the retrievals with abundant (2. Ox 10 -3 ) 
and scarce CO2 (10~ 7 ), and height-dependant CO2 profile. The 
functional derivatives for CO2 show low sensitivity at 0.1-1 mbar, 
which implies that the abundance of CO2 in the upper atmosphere 
remains underconstrained by existing data. 

Thirdly, we have quantified the degeneracies between the at- 
mospheric properties using the cross-correlation functions. The 
correlations between molecular abundances and temperature ex- 
hibit substantial degeneracies, in particular, between temperature 
and H2O abundance at 300 mbar. We determine the temperature 
uncertainty at different pressure levels by calculating the statistics 
of the retrievals with the variation of each molecular mixing ratio. 
In the same way, the molecular abundances are also constrained by 
demonstrating that the uncertainty on retrieved parameters, based 
on the available measurements and the degeneracy with the tem- 
perature profile is considerably larger than previous studies have 
suggested. Therefore, additional data is clearly required to break 
the degeneracy between temperature and composition, particularly 
H2O. As long as the number of retrieval variables for the molec- 
ular abundance are kept to a small number, the cross-correlation 
functions are not significant between molecules. If we were forced 
to use a more detailed representation of the vertical distributions 
of the gases and scattering properties, then the degeneracies would 
become more substantial. 

This study has shown the benefit of considering a wide spec- 
tral range for breaking the degeneracy between different atmo- 
spheric parameters. By constraining the overall shape of the in- 
frared spectrum, the HST and Spitzer secondary eclipse data al- 
low us to evaluate the vertical temperature structure and to place 
crude estimates on the molecular abundances in this hot Jupiter. 
Studies that focus on smaller wavelength ranges or fewer data 
points are subject to broader uncertainties than are currently being 
presented in the literature. Future work will permit an extension 
of this technique into cross-comparison with transmission spec- 
troscopy, including the visible range (alkali metals, metal oxides, 
and clouds and hazes) especially at the terminator regions of this 
planet. 
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